STABILITY OF A FORCE-BASED HYBRID METHOD IN THREE 
DIMENSION WITH SHARP INTERFACE 



JIANFENG LU AND PINGBING MING 

Abstract. We study a force-based hybrid method that couples atomistic 
model with Cauchy-Born elasticity model with sharp transition interface. We 
identify stability conditions that guarantee the convergence of the hybrid 
scheme to the solution of the atomistic model with second order accuracy, 
as the ratio between lattice parameter and the characteristic length scale of 
the deformation tends to zero. Convergence is established for the three dimen- 
sional system without defects, with general finite range atomistic potential 
and simple lattice structure. The key ingredient of the proof is regularity 
and stability analysis of elliptic finite difference schemes. We apply the result 
to an atomistic-to-continuum scheme for a triangular lattice with harmonic 
interactions. 



1. Introduction 

Multiscale methods couple together atomistic and continuum models have re- 
ceived intense investigations in recent years (see e.g. [TJ[TUl[2TJ[2Hl|3n]). Generally 
speaking, there are two main categories of methods coupling atomistic and contin- 
uum models: energy-based and force-based methods. The energy-based methods 
employ an energy that is a mixture of atomistic energy and continuum elastic 
energy. The energy functional is then minimized subject to suitable boundary con- 
ditions to obtain the deformed state of the system. The force-based methods work 
instead at the level of force balance equations. The forces derived from atomistic 
and continuum models are coupled together. The force balance equations supple- 
mented with suitable boundary conditions are solved to obtain the deformed state 
of the system. 

From a numerical analysis point of view, the key issue for these multiscale meth- 
ods is the consistency and stability analysis of the coupled schemes [lOl Chapter 7] . 
In this paper, we study force-based atomistic-to-continuum hybrid methods in two 
and three dimension with sharp transition between the atomistic and continuum 
regions. In our previous work [5D], we developed the stability analysis in general 
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dimension for a force-based atomistic-to-continuum method with smooth transition 
between the two regions. The main focus of the current paper is to extend the sta- 
bihty analysis of |20j to hybrid schemes with sharp interface between atomistic and 
continuum models. 

Comprehensive reviews for force-based hybrid methods can be found in j23} 
Section 5 and Section 6] and [29l Section 12.5]. A class of force-based methods 
uses a handshake region (transition region) between the atomistic and continuum 
regions. A representative of such methods is the concurrent AtC coupling method 
(AtC) developed in a series of papers [3ll4l [T3l[25] , blending the continuum stress and 
the atomistic forces. The method recently proposed and analyzed by the authors 
in [20] shares certain common traits with the AtC method. It is proved to be stable 
and convergent with optimal convergence rate. The numerical implementation of 
the method can be found in [33) . 

As a representative for force-based methods without handshake region, the FEAt 
method of Kohlhoff, Gumbsch and Fischmeister [17] is perhaps one of the earliest 
force-based methods. In this method, an elasticity model is coupled with an atom- 
istic model. The FEAt method does not use a handshake region: the transition 
between the atomistic model and the continuum model is sharp. This kind of cou- 
pling is generalized in CADD method 27' , which uses the discrete dislocation model 
in the continuum region. 

One of the main advantages of the force-based methods is that consistency is 
achieved with fairly simple construction. Hence the main focus of analysis of such 
methods is stability. The stability for one-dimensional force-based method was 
analyzed in [8l[9]. The generalization to higher dimension is quite nontrivial due 
to the complications of lattice structures and atomistic interaction potentials. The 
main idea in our previous paper |20| and the current paper is to establish linearized 
-ff ^-stability of the hybrid scheme by viewing the scheme as a nonlinear elliptic finite 
difference system and applying the elliptic regularity estimates for such system. 
The recent work [19j by Li, Luskin and Ortner proved linearized if^-stability for 
methods with smooth coupling under certain stability conditions. These conditions 
however were not yet known how to check explicitly [19]. They also studied how 
the size of the transition region affects the stability. 

For atomistic-to-continuum hybrid method with sharp interface studied in this 
paper, stability might fail at the interface. Therefore, to make sure that the hybrid 
scheme is convergent, we need to check the stability conditions at the interface 
for the coupling schemes. We shall identify the interface stability conditions as 
analogs of the complementing boundary conditions in elliptic PDE literature. From 
a physical perspective, these stability conditions amount to check whether there 
exists nontrivial surface phonon at the interface of the two schemes. To some 
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extent, these stabihty conditions are analogous to the famous Gustafsson-Kreiss- 
Sundstrom stabihty conditions )14] for finite difference approximations of mixed 
initial/boundary value problems. 

The main result in this paper is the linearized iJ^-stability and convergence of 
atomistic-to-continuum hybrid method under the stability conditions. The essential 
ingredients are regularity and stability analysis of finite difference schemes. As a 
consequence of our main results, we will show that a force-based atomistic-to- 
continuum coupling for a 2D triangular lattice with next-nearest neighbor harmonic 
interaction is stable and hence convergent, when the interface between the atomistic 
and continuum regions is parallel to the (l,-\/3)/2 direction of the lattice. Let 
us finally remark that while we focus on hybrid schemes coupling atomistic and 
nonlinear elasticity models, the ideas and techniques in the current paper can be 
extended to other force-based hybrid methods, e.g., the force-based coupling of 
peridynamics and nonlinear elasticity proposed in j26l . 

Before we formulate the problem and state our main result, we introduce some 
preliminaries and notations. 



1.1. Lattice function and norms. We will consider only Bravais lattices in this 
work, which is denoted as L. Let {aj}j^^ C M'* be the basis vectors of L, and d be 
the dimension, 

j 

Let {bj}'^^i C R'^ be the reciprocal basis vectors satisfying aj • bk — 2nSjk, where 
Sjk is the standard Kronecker delta symbol. The reciprocal lattice L* is then 

L* = |a; G M"* x = rijbj, n e Z^'j. 
j 

We take a computational domain 



and let ilg be a grid mesh in fl with mesh size e — 1/{2N), N £ Z 



+ ■ 



Using the reciprocal basis {bj}, we define 



1^ = ^ kjbj k e U^, -N < kj < ivj . 



We will identify functions defined on fl^ with their periodic extensions in this work, 
i.e., we consider the periodic boundary condition. General boundary conditions 
will be left for future work. 
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For /i e Z"^, we define the translation operator T^' as 

{T^u){x) = u{x + efijaj) for x e R'^, 

where the index summation convention is used. We define the forward and backward 
difference operators as 

= e-\T,^ - T) and D'^ ^ e'^I - Tf^^), 

where / denotes the identity operator. We say a is a multi-index, if a G Z'' and 
a > 0. We will use the notation |a| — '^j- For a multi-index a, the difference 

operator D" is given by 

d 

i=i 

where are the canonical basis of M.'^ (columns of a d x d identity matrix). 

We will use various norms for functions defined on fJ^. For integer fc > 0, define 
the difference norm 

\\<k- E J2 \{D^u){x)\\ 
0<\a\<k xen^ 

It is clear that \\-\\^ ). is a discrete analog of Sobolev norm associated with _ff'^(f2). 
Hence, we denote the corresponding spaces of lattice functions as H^{il) and L'^{^) 
when k = 0. We also need the uniform norms on il^, which is given by 

= max|u(a;)| , 
\\u\\^..^= E max|(i^»(x)|. 

0<|ct|<fc 

Recall that we identify lattice function u with its periodic extension to function 
defined on eL, and hence differences of the lattice functions are well-defined. These 
norms may be extended to vector-valued functions as usual. For k > d/2, we have 
the discrete Sobolev inequality IjitH^oo ^ ll^lle k- Here and throughout this paper, 
we denote A < B ii A < CB with C an absolute constant. 

The discrete Fourier transform for a lattice function u is given for ^ G L* by 

(1.1) ^2(^)^(^)'^e-«-„(:r). 

By the Fourier inversion formula, for a; € $7^, 

(1.2) u{x) = E ^"■^"(^)- 

We will use a symbol introduced by Nirenberg in [24j , which plays the same role 
for the difference operators as A^(^) = 1 + Aq(^) = 1 + |^|^ for the differential 
operators. For e > 0, let 
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and 

Auo = 1 + AL(o = 1 + E ^uo = 1 + E i s"^' 

It is not hard to check for any f e L* , there holds 

(1.3) cA^o < A?(e) < A2(e), 

where the positive constant c depends on {bj}. 



1.2. Atomistic model and Cauchy-Born rule. We consider classical empirical 
potentials: For atoms located at {j/i, • • • tUk}, the interaction potential energy 
between the atoms is given by V{yi, • • • , ijk), where V often takes the form 

Viyir-- ,yK) =E^2(2/,:/£,2/i/e) + E^3(y»/£,2/i/e,yfc/e) + -- - , 

where we have omitted interactions of more than three atoms. 

As in our previous work [20] . we will make the following assumptions on the 
potential function V: 

(1) 1^ is translation invariant. 

(2) V is invariant with respect to rigid body motion. 

(3) V is smooth in a neighborhood of the equilibrium state. 

(4) V has finite range and consequently we will consider only interactions that 
involve a finite number of atoms. 

For simplicity of notation and clarity of presentation, our presentation will be lim- 
ited to potentials that contain only two-body and three-body potentials, and we 
will only make explicit the three-body terms in the expressions for the potential. 
It is straightforward to extend our results to potentials with interactions of more 
atoms that satisfy the above assumptions by following the discussion on the three- 
body terms. As the potential function y is a function of atom distances and angles 
by invariance with respect to rigid body motion, we may write 

V3{yi,yj,yk) = V3 (^\y, - yjf' , \yi - yuf' , (y» - yj,yi - yk)^ , 

where (•, •) denotes the inner product over W^. 

We assume that the atoms located at VL^ are in equilibrium, with a; G f2e the 
equilibrium position. Positions of the atoms under deformation will be viewed as a 
function defined over fi^, which is denoted as y{x) = a; + u{x). Hence, m : fig — > 
is the displacement of the atoms. Define the space of the displacement field as 

Ae = < M : Sle ^ K'' I E "(^) = ^ f ■ 
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The atomistic problem is formulated as follows. Given force field fe : 
find u e such that 

(1.4) u = arg min /at(u), 

uex^ 



where 



and 



Here 5 is the set of all possible (si,S2) within the range of the potential. By our 
assumptions, S" is a finite set. We only make explicit the three-body terms in the 
potential as mentioned before. In /at, e'^ is a normalization factor, so that /at is 
actually the energy of the system per atom. 

The Euler-Lagrange equation for the atomistic problem is 

(1.5) T^t[u]{x) ^ feix) xen^, 

where 

J'.tMix) = E {^s, (29iy(,,^,,)[y](x)i?+y(x) + d3V^s,,sM{x)DXy{^)) 
(si,s2)es 



+ D-^ (292y(,,,,,)[y](x)i?iy(a.) + 93^^(s„s.)M(2;)/?+?/(x))j, 
where for i = 1, 2, 3, we denote 

the partial derivative with respect to the z-th argument of V . 

To guarantee the solvability of (|1.5p . we assume that takes the following form: 

/,(x) = £-'*/ f{z)dz, xel^e, 

Jx+eF 

where f{x) is a function defined on 51 with zero mean. This makes sure that /e(x) 
satisfies 



E fe{x)^e-^ f f{x)dx = 0. 



To introduce the Cauchy-Born elasticity problem [51[TT1[T2], we fix more nota- 
tions. For any positive integer k, we denote by W'''P{Q;R'^) the Sobolev space of 
mappings u: — )• R'' such that ||u||^yfc,p is finite, and by W^^ the Sobolev 
space of periodic functions whose distributional derivatives of order less than k are 
in LP{n). For any p > d and to > 0, we define X as 



X = iu-.n 



u 
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The Cauchy-Born elasticity problem is formulated as follows. Find u & X such 
that 

(1.6) M = arg min /(m), 

where the total energy functional / is given by 

I{u)^ [ (WcBiyuix)) - f{x)uix)) dx. 
Jn 

Here the Cauchy-Born stored energy density Wcb is given by 

' (si,s2)es 

where for AeR'^'"', 

W(s,^,,){A) = V (|si + siA\^ , \s2 + S2A\^ , (si + siA, S2 + saA)) . 

The range S is the same as that in the atomistic potential function. 
The Euler-Lagrange equation for the Cauchy-Born elasticity model is 

(1.7) Tcb[u]{x) = fix), 
where 

Tcb[u]{x) = -div(i:>AVl^cB(Vw(x))). 

Here DaWcb{A) denotes the derivative of VFcb(^) with respect to A. 

Since we are primarily interested in the coupling between the atomistic and 
continuum models, we will take the finite element mesh 7^ as a triangulation of 
rig with each atom site as an element vertex. The triangulation is translational 
invariant. The approximation space is defined as 

X,^{ue W^Pin-M") I u\t e Pi(t), VT e r,}, 

where Pi{T) is the space of linear functions on the element T. We denote by J> the 
force from finite element approximation of Cauchy-Born elasticity problem (|1.6p . 

1.3. Formulation of force-based hybrid method with sharp interface. We 

are ready to formulate the force-based hybrid method. We take a continuum region 

= l^ajjOj < xi < 1/2, < < 1, j = 2, - • • 

and denote g the characteristic function associated with ^l^'. g{x) = 1 if x e flc- 
fla ~ fl\flc ~ \x I q(x) = 0} is the atomistic region. Note that the continuum 
and atomistic regions are separated by two hyperplanes \x\ = 0} and \x\ ~ 1/2} 
as a result of periodic boundary condition. The simple geometry here is chosen 
for simplicity of presentation, using localization techniques as in |18| (see also the 
proof of Theorem 12] below), the analysis can be generalized to any J7c with smooth 
boundary. 
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We consider a force field defined as 

(1.8) ^hyM(2:) = (1 - 9{x))J'i,t[u]{x) + g{x)J',[u]{x), x G ft,. 

Due to the choice of g, in the atomistic region ila, the force acting on the atom 
is just that of atomistic model, while in the continuum region flc, the force is 
calculated from finite element approximation of the Cauchy-Born elasticity. Since 
g is taken to be the characteristic function, we consider here a hybrid method with 
sharp interface, i.e., there is no transition or buffer region between the atomistic 
and continuum regions. 

Given a loading /e, we find u G such that 

(1.9) {n,T,,y[u])ix) = f,{x) xen„ 

where for a lattice function g, llg projects g to a function with zero mean: 
{W,g){x): ^ g{x) - 9{x'). 

x' 

As in [30], the convergence of the hybrid scheme is tightly connected with the 
linear stability of the finite difference scheme. Thus, we will study the linearized 
operator of J^hy Let us denote "HhyM the linearization of J-hy at state u: Hi^ylu] = 
-j^, so that "Hhyiw] is a linear operator acting on a lattice functions w, which is 
given by 

Hhvlulw = lim \u + tw]. 
t-i>o dt 

We will rewrite the operator Hhy in the form of a difference operator as 

'^hyN = ^ hi,y[u]{x,fl)T'^, 

where the coefficient /ihy[i*](a;, /i) is a c? by c? matrix (probably asymmetric) for each 
X and fi Cz A, which is given by 

^{T^^u)|3{x) 

where a, P — 1^ ■ ■ ■ ,d are indices. Here A is the stencil of the difference operator, 
which is finite by assumptions on the atomistic potential. By the definition of Jiiy, 
we have 

(1.11) /ihyM(a;,^) = [l - g{x))h^t[u]{x,^i) + g{x)h^[u]{x,^i), 

where ft.at[u] and h^[u\ are given by similar equations as (jl.lOl) with J^hy replaced 
by Fat and respectively. 

Define h\iy[u]{x, ^) as the symbol of the pseudo-difference operator Hhy M, which 
is given by 

Ky[u\{x,Cj = ^ /ihyM(a;,Ai)exp(^«e^^jaj • for i e L*, 

tJ-^A j 
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and similarly for he [u] and h^t [u] . By definition, we have for any x G f^e , 

for j,k = l,...,d and similarly for he[u] and ft-atM- It is also clear that 
implies 

(1.12) hi,y[u]{x,0 = (1 - g{x))Kt[u]{x,0 + g{x)h,[u]{x,0- 

When we linearize J-hy around the equilibrium state u = 0, we will simplify 
the notation as Hhy — "HhyiO], ft-hy = ^hy[0], and similarly for those defined for 
atomistic model and finite element discretization of the Cauchy-Born elasticity. By 
the translation invariance of the total energy /at at the state u = 0, we observe 
that the coefficients of the symbols hi^tix,fj,) and he{x,fi) are independent of the 
position X, i.e., 

We also denote ?^cb as the linearization of Jx;b at the equilibrium state u = 0, 
and define /icB(a^,^) as its symbol. Due to the periodic boundary condition, the 
argument ^ in the symbol hc^{x,S^) only takes value in L*. Again, due to the 
translation invariance of the total energy, the symbol /icb is also independent of x. 

An elementary calculation shows that the matrices hat{^), he{S,) and hence 
hhy{x,£,) are Hermitian for any e > 0. As in [25, we make the following stability 
assumption about the atomistic potential: 

Assumption A. The matrix ft.at(0 is positive definite and there exists a positive 
constant Oat such that for any e > and any ^ S L* , 

det/iat(e) >aatA2;^,(e). 

We also recall from [20l Lemma 3.3] that, as a consequence of Assumption \^ 
for sufficiently small e, the finite element approximation is also linearly stable. 

Lemma 1.1. There exist constants a, Eq > that for any e € (0,£o); the symbol 
/le(^) is positive definite and for any ^ G L* , 

dethe{0>aAl'^M. 

1.4. Main result. The main focus of the current paper is to establish convergence 
for the hybrid method with sharp interface as e — > 0. In [20], convergence was 
proved for any short-range interaction potentials when g is a smooth function, or in 
other words, when the transition region between atomistic and continuum regions 
is of 0{e~^). In this paper, as the transition is sharp, we require additional stability 
conditions to guarantee convergence. Our main result is 
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Theorem 1 (Convergence). Under Assumptions\^\B^ andQ there exist positive 



Remark. The Assumption [XI is a natural stability condition for the atomistic lattice 
system. Compared with the convergence result in [20], the additional assumptions iBl 
andOcome from the coupling between the atomistic potential and the finite element 
discretization at the interface. The two assumptions will be given in Section [S] 

Remark. While we do not attempt in this work to optimize the regularity assump- 
tion on /, we note that it is easy to relax the assumption to f £ W^'P{il.) with 
p> d. 

The proof of Theorem [1] follows a similar strategy as in ^20^ . Actually, once we 
obtain the stability estimate Theorem 21 the proof of Theorem [T] is essentially the 
same as that of [5D1 Theorem 1.1]. The consistency analysis of the scheme follows 
from that of |201 Section 2] with some immediate adaptations. Observe in particular 
that the proof of consistency does not depend on the smoothness of g. Hence, we 
will focus on the linear stability analysis, and omit the consistency part. 

The remaining of the paper is organized as follows. In Section [51 we recall 
the concept of pseudo-difference operators, which will be used to establish the 
regularity estimates in Section |3l The stability estimate then follows from the 
regularity estimate combined with consistency, which is presented in Section |4l In 
Section we apply the general theory to an example of force-based atomistic- 
to-continuuni hybrid method for 2D triangular lattice with next-nearest neighbor 
harmonic interaction. We make some conclusive remarks in Section [6l 

2. Regularity of finite difference scheme up to the boundary 

To analyze the general hybrid scheme, we will take the viewpoint of our previous 
work [5D] and regard the scheme as a nonlinear finite difference scheme. One of 
the main ingredients we will use in this paper is the regularity estimates up to the 
boundary for elliptic difference systems established in [28]. We also note that the 
regularity estimates of elliptic difference equations and systems have been investi- 
gated by several works [6l[T5 l [T6 l [T8 l l22 | l28 | l31j . For reader's convenience, we recall 
here the setup and results with some adaptation to the current work. 

2.1. Difference operators. We write a difference operator in the form of 




(2.1) 
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where the stencil AdU^ \s finite. We define the symbol of 

(2.2) U(xX)^Y.^e{x,y.)e^^'>^< 

for any ( £ 'E'^. By definition, we have 

We choose A as the minimal stencil, that is, for ii £ A, the symbol li;{x,fi) is not 
identically zero. For this choice of A, we let 

/i. = min /ii, /Ij — max /i^, i = 1, • • • , d, 

and call (/i.,/!^) the extent of the difference operator in direction. 

2.2. Elliptic difference system. Consider a half-space 

j 

with periodic boundary condition in y. i^haif is discretized by grid points {x^,y^) 
with mesh size e = 1/(2N). Here x^, = evai with v >Q and Vfj, = £ J2j i^j'^j+i with 

Consider a linear system of difference equations 

n 

(2.3) ^{Le,tjUj){xu,yt_i) = fi{xy, y^), i = 1, • • • , n 
with boundary conditions 

n 

(2.4) ^(Be,fejWj)(xo,y^) = 0fc(y^), = 

We denote l^^ij the symbol of L^^ij and the symbol of B^ ^j viewed as differ- 
ence operators on the whole space, defined as in (12.21) . For our purpose, it suffices 
to consider the case that the coefficients of difference operators are independent of 
X, y, and e, i.e., 

le,tj {x,y,iy,fi) = kj {v,fJ.), be,kj {y, V, fi) = bkj (j^, /i) . 

We will henceforth make this simplification. Note that the symbols le^ij and f'e,fcj 
however still depend on e. 

Definition 2.1 (Regular elliptic difference system). We call (|2.3p a regular elliptic 
of order (cr, t) for cr, t G Z" if the following conditions are satisfied, 
i) For each i, j = 1, ■ ■ ■ ,n and e sufficiently small. 



^To be consistent with 1201 . our scaling of the reciprocal space in terms of e is slightly different 
from that of 128. . 
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ii) There exist positive constants Co and eo such that 

|detZe(C)|>A,(C)2p 
for all < e < £o and \(\ > Co, where 2p — J2ii'^i + '''»)■ 

As a convention, we will choose a and r so that niaxdi — 0. 
We also take p e so that 

(2.5) |5,,,,(C)|<Ae(C^+^^ 

and pk is the smallest possible integer satisfying (|2.5|) . (p, t) gives the order of 
viewing as difference operators on the whole space. 

We will assume that the operators L^^ij and B^^kj only contain differences of 
order ct^ + Tj and pk + Tj respectively. 

Assumption 1. Let the difference operators Le,y have extent {j^^j , Vij ) in x— direction, 
we assume that there are a"*" , fi^ in Z'' such that 

< hLij < < a+ + i,j = !,■■■ ,n, 

and such that the number of roots z(?7), counting multiplicity, of the equation 
Riz,7j) = is + /?+). Here 

i?(z,77) = det ^ ^ hj{iy,p)z''e"'-'^ 

= det {-ie'^ log z, e"^?y)] . 

Furthermore, the number of boundary conditions q is equal to the number of roots 
of R{z, 77) = which satisfy 

< |z(?7)| < 1 for ^ 0. 

In the case where q is larger than p, we need an additional assumption. We 
assume that the boundary conditions are ordered so that pk are in increasing order. 
We define 

p= max (pfc + 1,0), 

l<k<p 

and 

minp/s + l, ifg>p; 

k>p 

CO, if g = p . 

Assumption 2. We assume that the boundary operators satisfy pk > p for p < 
k<q. 
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In the limit £ — )■ 0, we obtain from the system of differential equations associated 
with the system of difference equations. 

n 

^{L,j{d^,dy)uj){x,y) = fi{x,y) i = !,••• ,n, 

i=i 

n 

^{Bkjid^,dy)uj){0,y) = (t)k{y) k=l,--- ,p. 

i=i 

The symbol of Ljj , denoted as lij (C) , is given by 

We will make the following assumption on the symbol hjiC)- 

Assumption 3 (Supplementary condition). The determinant det[lij{()] is of even 
degree 2p. For every pair of vectors G K"^, the polynomial det[lij{( + t(')] in 
the complex variable r has exactly p roots with positive imaginary part. 

Let G be the tangential grid given by 

G = {y^\0<tij<2N,j = l,--- ,d-l}. 

We define discrete Fourier transform for lattice functions u with respect to y as 

for T] e Z'^-i with \T]j\< N fov j = !,■■■ ,d- 1. 
We define the following norms 



»)6L* ^ iy=0 m=0 



and 

i|2 



U|i:^.. = Ell"jlls+r„e 



s+T,e 



At the boundary, we define 



[si 

\^j\l,= E E|Ae('?)^"™^"e%■(^0,ry)|^ 



'?eL;,j m=0 

and 



I'^l'+r.e = El"J■|'+r„e• 
j=l 
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2.3. Complementing boundary condition. Due to periodicity in y, we may 
reduce tlie system to a one dimensional operator by Fourier transform. Recall that 

Define coefficients l^^ij for 7] e Z'^-i as 
The reduced operator is given by 

Analogously, we define B^^i^j as the reduced operator of Bs^kj ■ We denote B]. the 
operators corresponding to the first p boundary conditions and B^ for the remaining 
q — p ones. 

Similarly, we define the reduced operator for differential operators as 

Lij{dx,6) = Lij{dx,ie) 

for 9 & S'^, and analogously for Bkj. 

For the reduced operators, we consider three types of eigensolutions, which are 
defined as follows. 

Definition 2.2. We call an eigensolution of type I a nontrivial solution of 
{Ls{ri)w){xv) =Q for some 77 7^ 0, 

satisfying 

(a) (Bi(r7)u;)(xo)-0; 

(b) (Bf(r7)w;)(xo) = 0; 

(c) wi^Xv) as 1/ ^ DO. 

Definition 2.3. We call an eigensolution of type II a nontrivial solution of 
{L{d^, e)w){x) = 0, for some 6 e S"^, 

satisfying 

(a) {B\d^,e)w){x^)^Q- 

(b) wi^x) — >■ as a; — > 00. 

Definition 2.4. We call an eigensolution of type III a nontrivial solution to 

{LMw){x) = 0, 
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satisfying 

(a) {B^,{0)w){xo) = 0; 

(b) ^(a;^) — ^ as — > oo. 

Remark. We remark that type I and type III eigensolutions are connected to the 
notion of surface phonon in the physics hterature. The complementing conditions 
give a mathematical characterization of the surface phonon. On the other hand, 
the type II eigensolutions are surface waves at the boundary of the PDE. 



Definition 2.5 (Complementing Boundary Condition). The system (|2.3|) 
with boundary conditions ()2.4[) satisfies the Complementing Boundary Condition if 
there are no eigensolutions of type I, II, or III. 

The main result of [28l Theorem 3.1] is the following regularity estimate of the 
solution of elliptic system with boundary conditions. 



Theorem 2. If u is a solution to the system (j2.3p with boundary conditions (12. 4p 
and Assumptions [3 and O are satisfied, then the following regularity estimate 
holds for each s with p < s < p* and e sufficiently small, if and only if, the 
Complementing Boundary Condition holds. 



\Ml+s.e + \^\l+s-l/2.e < \^i\l-p-l/2,e 

(2.6) 



where t = p + ^ [2{s — p)\ . 



+ \\f\\l-a.s + \ML^ 



3. Regularity estimate for hybrid schemes 

In this section we will apply the regularity estimate Theorem [5] to atomistic-to- 
continuum hybrid schemes. 
We consider domain 

f^strip = [xai +J2yj"-j+i I a; e M, 2/ e [0, l)"*"^! 

j 

with periodic boundary condition in y variable, i^strip is discretized by grid points 
{xij,y^) with Xi, = evai with v € 1, and — £j2j l^j'^j+i with p, e Z''^^. We 
consider the following hybrid system on fistrip: 

d 

(3-1) ^{'He,ijUj){x,,,yf,) = f,{x^,yf,), i/ < 0, z = 1, • • • , d; 

d 

(3-2) ^{niit,ijUj){x^,yf,) = fi{x^,yf,), > 0, i = 1, • • • , d. 

j=i 

We will reformulate the hybrid scheme as a system of difference equations with 
boundary conditions by folding with respect to the surface {x — 0} (similar folding 
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trick was used in [7]). Let and (Eij^^fj) be the extent oi'He.ij and 'Hat,ij 

in X direction respectively. To simpUfy the presentation, we wih assume that the 
extents are the same for aU i,j's, which is usually the case for applications in 
atomistic-continuum hybrid schemes. We denote them as (i^'^,!?'^) and {v^.v"")- 
The construction can be extended to more general cases, which will be omitted for 
simplicity. Without loss of generality, we also assume that i/^ < i/^ < Q <V'^ <V°^ . 
Let us rename the variables as 

Ud+i{xy,yfj) = Ui{xjja+^,y^), 1^ >0,iie W^'^. 

This leads to d{T^ — i/^) compatibility conditions 

(3.3) C/i(a;^,?/^) = C/d+j(a;i7<=-^a-,.-i, y^), < <V'' - j/" - 1, i ^ 1, ■ ■ ■ ,d. 
We will rewrite equations p.ip - p.2p in terms of U's. Define 



Le,ij 



X] KtXi-d)(3-d){v,n)Tl^T^y, i,j ^ d+1,- ■ ■ ,2d; 

0, otherwise. 
By construction, the extents of the operators are given by 

iE,j,T^v) = I (0,17'^ - z."), z,i - d + 1, • • • , 2d; 
[(0,0), otherwise. 

Let 

/i(a;_^_i,y^), i = !,■■■ ,d; 



fi-dix^,y^,), i = d+l,--- ,2d, 



we then have 

2d 

(3.4) ^(L,,,jC/,)(x„j/^) = F,(x„2/^), iy>0,^ieZ''-\ 
i=i 

We further define for fc = 1, • • • , d{T''^ — v^) and j = 1, • • • , d, 

where? = [(fe— l)/(ij and I = [{k—1) mod d] + l. The compatibility conditions 
are then equivalent to the boundary conditions 

2d 

(3.5) ^(B,,fe,C7,)(a:o,y^)=0, ^i eZ''-\k = I,- ■ ■ ,d{v- - v^). 

Therefore, we have reformulated p.ip - p.2l) into a difference system p.4p with 
boundary conditions p.Sp . 
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Since 

(3.6) |det;,,,j(C)| - |det/i,,,j(C)| \det h,t,rj{C2)\ ■ 

The regular ellipticity of the operator is guaranteed by Assumption El as Hat 
and He are both regular elliptic by Lemma [TTTl Moreover, the order of is (2, 0), 
and hence, p = 2d. 

Taking a+ = 0, /3+ 17^ - i/^ j = 1, • • • , d, and /3+ = 17'^ -z^", j = rf+ 1, • • • , 2d, 
we have then < v_^j < Vij < af + . The number of boundary conditions we 
have imposed on (I3.5P is g = d{V'' — v^)- Define 

R'^{z,i]) = det [he^ij{—ie~^ log z,e~^ri)] ; 
R"^{z,r]) =det [hi,t,ij{-ie'^ log z,e~'^ri)] . 

From (I3.6p . the Assumption [T] for the system (|3.4p - (|3.5l) is equivalent to 

Assumption B. The number of boundary conditions (i(l7^ — e'^) is equal to the 
total number of roots of R^iz, r/) — and R"'{z, rj) = which satisfy 

< |z(r/)| < 1 for ^ 0. 

We remark that this assumption is satisfied for linearized finite difference oper- 
ators given by common inter-atomic potentials. 
By definition, we have 

Pk = [{k - l)/d\ - 2. 

Hence, p = and p* — 1. Assumption [2] is satisfied by the system p.4p - (|3.5p . 

The Supplementary Assumption |3] is automatically satisfied since 2d > 3 as a 
result in [2]. 

The stability of the hybrid scheme relies on the stability at interface, which is 
characterized by the complementing boundary condition. 

Assumption C. The hybrid system (|3.ip - p.2p . or equivalently (|3.4l) - (|3.5p satisfies 
the Complementing Boundary Condition, i.e., there are no eigensolutions of type 
I, II, or III. 

This stability assumption needs to be check for particular atomic interaction 
potentials. We will discuss an example in Section [5l 

We can now state the regularity estimate for the hybrid system p.H) - (|3.2p . 

Proposition 3.1 (Regularity). If u is a solution of the system (|3.ip - (|3.2p . under 
Assumptions\^\B^ and[3 the following regularity estimate holds for e sufficiently 
small 



(3.7) 



ll"ll2,.^ll/llo.e + ll"ll 
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Proof. Applying Theorem [2] to the system p.4p with boundary conditions p.5|) . we 
have 

ni2,e<ll^llo,. + nio.e- 

By definition of F and U, it is clear that \\F\\q^^ = ll/llo,6' ll^^llce ^ 2 and 
\\u\\2 e < 11^112 £• Therefore, we obtain 

\\u\ke<\\f\\o,e + \Mo,e- 

This completes the proof. □ 

The regularity estimate for the linearized operator Hhy = "HhyiO] then follows. 

Theorem 3 (Regularity). Under Assumptions\^ [21 a^irfLO for any v E H^{^1), 
we have 

(3.8) IIHl2,e^ll^hyf|lo,e + l|f|lo,e 

for e sufficiently small. 

Proof. The theorem follows from Proposition 13.11 combined with a standard local- 
ization argument as in |18| . We sketch the proof here for completeness. 

Let Lpi and ip2 be a partition of unity of the domain 17, such that if if 2 € C°°(r2), 
+ (^2 = 1, ^i{x) = 1 if e [3/8,5/8), and ifi{x) = if xi e [0,1/8) U [7/8,1). 
We then have 

M\2,^ < Il'^ll'll2,e + ll'/'2l'|l2,e- 

The estimates for the two terms on the right hand side are the same, we only 
consider [[(piiiHj^ in the following. Let u = fiv, since ipi is compactly supported 
in ri, we can extend u by as a function on the grid points in r^strip- Note that 

ll'^hy'"llo,e = Il'^hy(<y3iw)|| < HHhyfjla^^ + l|t'llo,e > 

since ipi is a smooth function. Thus, by Proposition 13. 1[ we have 

\\'PM\2..e - ll'^hyU||o,£ + \\u\\o,e - Il'^hy^'ll2,e + ll^llo,e ■ 

Therefore, we obtain the regularity estimate p.8p . □ 

4. Stability 

The main theorem in this section is the following stability estimate. 

Theorem 4 (Stability). Under Assumvtions \M\^ and\^ for any v E H^{^) with 
mean zero, we have 

(4.1) IK'IU2 ^ ll^hyl'll.^o 

for e sufficiently small. 
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To obtain the stability estimate from the regularity estimate Theorem [3l we 
need to eliminate HfH^Q on the right hand side of p.Sp . The proof is based on 
the uniqueness of the continuous system from ellipticity, the consistency of the 
finite difference schemes to the continuous system, and the regularity estimate 
Theorem [3l We remark that the stability analysis in our previous work [20j relies 
on the smoothness of g, so it does not apply to the sharp transition case. The 
current proof does not require smoothness of g and so applies both to the hybrid 
method with smooth or sharp transitions. 

In order to connect the finite difference system with continuous PDE, we need to 
extend grid functions on fig to continuous functions defined in fl. For this purpose, 
let us define an interpolation operator as follows. For any lattice function u on 
r^e, we define Q^u G L^(ri) as 

(4.2) iQeu){x) = ^ e O. 

Comparing with (jl.2l) . we know that Q^u agrees with u on fJ^. We have the 
following properties of Qe ■ 

Lemma 4.1. For k > 0, there exists constants Ck, Ck > 0, such that for any u, 

Ck I|w||^fc(f2) < IIQeUjl^fc^j^j < Ck ■ 

Proof. It is not hard to check for any ^ e L* , it holds 

cA2(C) < A2(C) < K\0- 

The conclusion follows. □ 

Let X be a standard nonnegative cut-off function on W^, which is smooth and 
compactly supported, with = 1- Let Xe be the scaled version 

Xe(2:)=e-("^'x(£-":r), 

for some a with < a < 1. The choice of the value of a will be specified in the 
proof of Proposition 14.41 

Define a low-pass filter operator K,, for / e L^(il) using Xe as Fourier multiplier: 

K7f{i) = (2^)^/(0S(0 = (2^)'^/(0x(£"C). 
In real space, convolves / with Xe- Integrating by parts, we obtain 

(4.3) \Xem<Ck\e"£.\-\ foraUfceZ+, 

(4.4) (27r)'^S(0) = 1. 

Hence, is indeed a low-pass filter. For simplicity of notation, we denote 

Ue = KsQeUe, 

for lattice function on fi^. 
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Let US first recall some consistency results proved in [201 Section 2] (proofs of 
these results do not depend on the smoothness of g) , which will be used in the proof 
of stability below. 

Lemma 4.2 (Consistency). For any u smooth, we have 

(4.5) ll^atM --^cbMIIl^ < Ce' \\u\\^^s.^ , 

(4.6) \\:F,[u] - ^cbMWl^ <Cs 



(4.7) II J-hyN - -FatNllioo < Ce' hlUis,oo , 

where the constant C depends on V and \\u\\j^ao, but is independent of s. 

Using Lemma l4.2[ we prove the consistency results for linearized operators. 

Lemma 4.3 (Consistency of symbols of linearized operators). There exists eo > 
and s > such that for any s < Eq, x € o,f^d ^ S L*, we have 

(4.8) \Kt[0-hcBm<Ce\l + ny\ 

(4.9) \h.t{0-KyM\<Ce\l + \i\y\ 

Proof. The proof for (|4.8p and (|4.9p are analogous, and hence we will only prove 
the latter. By definition, for 1 < j, A: < d, 

(/iat),/c(0=e-"-«(Hat(efe/e)),(a;), 
(V)./c(2:,0=e-"-«(Hhy(efe/5)),(x). 

where f^{x) — e'^ '» for x G il. Taking difference of the above two equations, we 
obtain the bound 

\Kt[0-Ky{x,0\<C sup ||Hat(efc/«)-Hhv(efc/«)|li..- 

l<k<d 

Note that by the definition of linearized operators Hat and "Hhy, we have 
'Hat(efc/e) - 'Hhy(efe/^) = ^fim ^(J"at[i(eA;/e)] - J'hy [^(efc/e)]) ■ 

Hence, 

1 

t^o+ t 

< llefc/dlvK-.o. < ||e./dl^. < e'il + |^|^)^/^ 

where s is chosen so that the Sobolev inequality ||/||vi/i8.~(n) — ^ll/llff=(0) holds 
for any / G 7?'*(n) (s depends on the dimension). Here, we have used Lemma 
in the first inequality, noticing that ||tefc/c||j^oo is uniformly bounded for ^ as t — >■ 0. 
This concludes the proof. □ 



||'Hat(efe/^) -'Hhy(efc/5)||^„e = lim - WJ^e^Mekf^)] - -T^hy [^(efe/^)^^,. 
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With these preparations, we now state a key proposition wih be used in the proof 
of Theorem [H 

Proposition 4.4. For {w£}e>o such that Ve G H^{^) and \\ve\\^2 uniformly 
bounded, we have 

Proof. By triangular inequality, 

ll'HcBWe - 'HhyWe||^2(o) - H'^hy'^e ~ "^hyWe || ^2(s;2) + H'^att^e - '^hyWe || ^2(q) 

+ II^CBWe - '^atW£||^2(o) • 

Hence, it suffices to show each term on the right-hand side goes to zero as e — 0. 
Step 1. By Lemma |4. II and the definition of 'Hhy, we have 

ll'Hhy^^e ~ '^hy1'e||^2(Q) < W^-hyVe " '^liy^^ell^.o 

We now show that WHeVe — T-L^VeW^ o converges to zero as e — > 0, the argument for 
the other term is identical. By Parseval's identity, 

(4.10) Wn.V, - Hel^elUo = WhemMi) ~ lR(Ollp(L;) ■ 

Note that \\he{C)Ve{C)\\i'i(i^*) = ll^e'^elleo bounded as e H^{^) and that 
Xe{£,) — > 1 for any ^ € L*, as e — 0. By dominance convergence, the right-hand 
side of (|4.10p converges to zero in the limit. Therefore, HHhyWe — ^hy'^e||2,2(Q) — ^ 0. 
Step 2. By Lemma 14.11 we have 

Hence it suffices to estimate the right-hand side. We calculate 
Therefore, using (|4.9p in Lemma [4.31 we have 

\\n.tv,-n^yVs\\lo<T.^'(^ + \i\'y\^M'- 

To estimate the right-hand side, note that by ()4.3p we have for any k E Z+, 

Taking k sufficiently large so that (1 + |^|^)'* ICP^*^ is bounded for ^ G L*\{0}, we 
then have 

\\H.tVe-m.yVe\\lo<s'e-'"^' 

Choosing a in the low-pass filter operator to satisfy 

(4.11) a < 2/k 
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guarantees that [["Hatt^e — '^^y^^W L^{n) ~^ ^is e — >■ 0. 
Step 3. By definition, 

n^eio = /icB(e)x(£"eR(e)- 

For the discrete system, we have 

Hence, we have, combined with (|4.8p and (|4.3p . for some s > and any fc G Z4 

<Cfc£4E(i + i^i)'i="^i"''i^'^-(^)i'- 



A—2ak |L, ||2 



We choose a sufficiently large A:, so that the right hand side is bounded by Ce 
The choice of a in (|4.1ip guarantees that ||Hcbi^£ — ^ati^£||^2(f2) — ?> as e — > 0. □ 

Proof of Theorem j^} Suppose (14. ip does not hold, then there is a sequence of func- 
tions {wk} and £a; > such that 

IIwajII^^ 2 ^ as fc — > 00; 

E Wk{x) = 0, for all A;. 

Set Vk ^ Wk/ \\wk\\^^ 2' '^'^ then have 

(4.12) ll«fclie.,2 = l for all fc; 

(4.13) ll'^hyWfell^^ 0, asfc^-oo; 

(4.14) E Mx) = 0, for all fc. 



Since 



Since HHhyUfcH^^ ^ 0, we have 

||'HhyWfc||^2(n) ^ 0' asfc-^cx). 
Moreover, by Proposition 14.41 

ll'HcBWfc - 'HhyWfc||^2(o) ^ 0' asfc^cx). 

Hence ||'HcBWA:|li2(-Q') — 0. Note also that the average of Vk is zero, since Vk{0) — 
0. By the invertibility of "Hcb on the subspace orthogonal to constant function. 
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Wfc||^2(Q) — >■ 0, as /c — >■ cx), while H^'feHej^ 2 — follows then ||wfe||^^ ~^ Indeed, 



since 



Vk\ 



for any (5 > 0, there exist ^ > and ki, such that for any k> ki, 

(4.15) < ^Z^- 

On the other hand, due to (|4.4I) . there exists ^2, such that for k > k2 

(4.16) ^ ||«'^fe(e)VI«'^fc(0r|<'^/4- 

eeL;^,|C|<s 

Moreover, as ||wfc||^2 —J' 0, there exists fcs, such that for k > k^, 

(4.17) Yl l^fc(Or<'5/4. 

€eiL.*,,l£l<H 

Combined ()4.15p - ()4.17p together, we have for k > max(fci, ^2, k^), 

Hence, limk^ao Ikfclle^ ~ '^^ From Theorem |3l this implies 

lim ll"fcllefc,2 = 0- 

The contradiction with the choice of Vk proves the Theorem. □ 

Theorem |4] extends to a deformed state it, we omit the proof which is the same 
as EQl Theorem 4.61. 



Theorem 5 (Stability). Under Assumptions\^\^ and\^ there exists 5 > Q, such 



that for any £ > and u, ||w||^^2,oo < 5 and any v € iJ^(f2) with mean zero, we 



have 

(4.18) \\v\l^^<C\\Ii,'H^M^\lfl^ 

where the constant depends on 5, but is independent of u, v and e. 

5. Examples 

We apply the result here to a particular example of atomistic-to-continuum hy- 
brid method. Consider a force-based method on a triangular lattice. Figure [T] gives 
the geometry and the choice of the atomistic and continuum regions. The interface 
between the atomistic and continuum regions is parallel to the (l,-\/3)/2 direction 
of the lattice. The interaction between atoms is assumed to be harmonic (quadratic 
potential). The force balance equation reads 
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Figure 1. Illustration of the geometry and coupling interface of the 
hybrid method. The coordinate of each atom is {v,f^), and the inter- 
face r along the atoms labeled with (0, /i). The interface F divides the 
domain into the atomistic region Sla and the continuum region Q,c- 



where J^hy is given by 

Ti,t[z]{x), X e ^a, 



with 



[ J^e[z\[X), X e Uc, 

1 

J"at[z](a;): = — ^[z{x + Si) - z{x)], 



i=l 
6 



Te[z]{x): = -^'^[z{x + s^) - z{x)], 



' ^=l 

where {si}}li is the interaction range of the first and the second neighborhood 
interaction of the triangular lattice. 

We rewrite the coupled force balance equation into a system. Define the map 

Tx = (-a;i,X2), 

and denote 

y{x) = y{Tx), and z{x) = (y(x), y{x)) . 
The hybrid method can be written as 

(5.1) J^hy[z]ix)^l,{x), 

where 

and e,{x) ^ iMx)J,iTx)). 

The equilibrium equations (|5.ip are supplemented with the following boundary 
condition 

(5.2) z{x) = Sz{Tx) on T, 
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where F is the interface as shown in Figure [Hand 

S = 



r 

1 0, 



Applying Fourier transform in the tangential direction, we obtain 
^e = ^ {Ts, (1 + e-'«-^^ ) + r_,, (1 + e'«-^= ) + 2 cos(C • S2) - 6) , 



+ T2s,e-'^-'' + T^2s,e'^"^ + 2cos(e • S2) - 12 

The first step is to consider the distribution of the root of the following two 
characterization equations. Let C, = e*^'*^ = e*^, obviously C 7^ 1 unless 9 = Q. 

(5.3) z(l + C) +2-1(1 + C) + C + C-6 = 0, 

and 

z'^C + z-^C + z{i + C + C + C'') 
^^'^^ +z-i(i + C + C + C') + C + C-i2 = o. 



The equation (|5.3p can be written as 



2 6-C-C. 



^^z + C = 0. 

1 + C 

The roots zi and z^ of the above equation satisfy 1^122 1 = I CI = 1- We cannot have 
ki| = 1-^21 = I7 otherwise det J> = 0, which contradicts with the bulk stability 
condition. Therefore, we have two distinct roots, one inside the unit disk, the other 
outside. 

To deal with ()5.4p . we let 2 = w^^/^, and write (|5.4p as 

(5.5) + + (Cl/2 + Cl/2 + (3/2 + C3/2)(^ + ^-1) + C + ^-1 _ 12 = 0. 

Let s = w + w^^ and A — ^^/^ _|_ ^1/2 ^ ^3/2 _|_ ^3/2^ write the above equation as 

+ As + C + C-U = 0. 
Denote by /(s) = s"^ + As + ( + C - 14:. Note that 
/(2) = 2A + C + C-10 

= 2(ci/2 + - 2) (^(c^/^ + C'/2 + 5/4)' + 

< 0. 
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It is clear that /(2) < unless C = 1, which is impossible. Proceeding along the 
same line, we obtain 

/(-2) = -2A + C + C-10 

= -2(Ci/2 + C'/2 + 2) + C'/2 - 5/4)' + ^ 

< 0. 

We immediately have, there exists two roots si and S2 with si > 2 and S2 < —2. 
This yields that there exists four roots {wi}j^^ of (|5.5p that satisfy wi > 1,W2 < 
1, — 1 < < and W4 < —1. Moreover, the roots satisfy 

W + W^^=Si, W = Wl,W2, 

W + W^"^ = S2, W — W3, Wi. 

To sum up, there exists four distinct roots for (j5.4p such that two inside the unit disk 
while the other two outside the unit circle. To sum up, we require three boundary 
conditions, which is consistent with (15. 2p . 

We next verify Assumption [Cl for the coupling scheme. 

For mode I, we have the following form of the solution 



z(x^,6')=ci z^ + cal Z2+C3LU3) \z^\<l, i = l,2,3. 




where zi denotes the root of (j5.3p with norm less than 1. Z2 and Z3 are the roots 
of (j5.4p with norm less than 1. As — >■ 0, we have \zi\ , \z2\ ~^ 1 while |z3| — )■ 3— 2-\/2. 
Substituting the above expression into the boundary condition, we obtain 



/ 1 


-1 


-1 


w 


— Z2 






The determinant of the matrix is 

{Zl - Z2){Z2 - Z'i){zs, - Zl) 



^0 



Z1Z2Z3 

because of the following lemma, whose proof is deferred to the end of this section. 

Lemma 5.1. zi,Z2 and Z3 are distinct roots. 

We conclude that there does not exist mode I eigenfunction. 
For mode II, notice that by definition, 

{HcBid^, 0)w) {x) = 6{dl - 0^)w{x) = 0, 

it is clear the only solution satisfying w{x) — > as x — >■ ±c» is the trivial solution 
w{x) = 0, hence mode II eigenfunction does not exist. 
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For mode III, we have C -> 1 as 6' — >■ 0. The solution takes the form 
z^c z\ with Z3 = 2\/2 - 3, 
substituting the above expression into the boundary condition, we obtain 

which yields c = 0. This concludes that there does not exist mode III eigenfunction. 
Therefore, the coupling scheme is stable and convergent. 

Proof of Lemma I5.il It is clear that 

Z2 = W2C^^^ and Z3 = WsC^^^ 

with — 1<W3<0<W2<1, this implies Z2 ^23- 
A direct calculation gives 



6-C-C-V(4-C-C)(2-C-C: 

Zl = = 

2(1 + 



6-c-c-v(4-c-c)(2-c-c) 

2(i + C)(i + C) ^ 



6-C-C- a/(4-C-C) (2-C-C) / N 

— V V ; / I 71/2 \ ^1/2 

2(2 + C + C) ^ ^ ■ 

Recalling ^ = e'^ with 9 e (— tt, tt), and we may write 

^ 2cos(g/2) ^,g/2 

^ 3 - cos 61 + ^(7-cos6')(l -cos6») 



Note that 

2cos(6l/2) 



3 — cos 6 + \/ {7 — cos d){l — cos ( 

this implies 21 7^ Z3. 

It remains to prove zi ^ Z2. Note that 



> > W3, 



with 



22 

2 



S = -^/2 + .^A2/4 + 14- (C + C)- 

Using 

A = C + C + C' + C'= (C + C) (C' + C') =4cos(0/2)cos0, 
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we write 

AV4+14- (C + C) = 16 cos2(6i/2)cos2 6* + 14 -2 cos 61 

= 16 008^(61/2) cos^ 61 + 14 - 2(2 cos2(6'/2) - 1) 
= 16-4cos2(6'/2)sin2 6l. 

This gives 



B^2^JA- cos2(6'/2) sin^ 6* - 2 cos(6l/2) 0086*. 
To prove zi ^ Z2, it remains to show \zi\ ^ |z2|, i.e., 



6 — COS 61+ cos6')(l — cos a 

Actually, we shall prove that for 9 E (— 7r,7r) and 0^0, there holds 



(5.6) , , 

2 V ^ 3 - cos 61 + ^(7- cose*) (1 - cos( 

The above inequality is equivalent to 

(5.7) 3 - cos 61 + v/(7- cos6')(l - cos 6*) > cos(6i/2) (b + \/B^ -A 
Denote by t = cos(6'/2), we write the above inequality as 

(5.8) 2-f^ + v/(4-<2)(i_<2) > i (^^(i) + ^g2(i) _ 1^ ^ i ^ [0, 1), 

where 

g(i): = t - 2t^ + 2v/l -t^ + f^. 
To prove (jS.Sp . we firstly prove 

(5.9) 2-t^>tg{t) te[0,l). 
A direct calculation gives 

2-f- tg{t) = 2(1 - + 2t (t^ - Vl - + t^) 

= 2(l-i2) + 



,2^ , 2t(t4-l) 



V]~t^Tt6 + t3 

2(1 -t^) ("i- 



t + i3 



VT^F+F + <3 



Note that 



VT^W+l^ > t, 

which follows from (1 — t'^){l — i*) > 0. Combining the above two inequalities, we 
obtain (lO)) . 

Next, by ()5.9p and note g{t) > 0, we obtain 

(4 - t^){l - f ) = (2 - -t^> t^{g'^{t) - 1). 

A direct calculation gives that g{t) > 1. Therefore, 



v/(4-t2)(l_t2)>i^^2(i)_l^ 
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which together with (|5.9p gives (|5.8p . This imphes zi ^ Z2 and completes the 
proof. □ 

6. Conclusion 

We have identified stabiUty conditions, especiahy stabihty conditions at the in- 
terface, for atomistic-to-continuum hybrid methods with sharp interface. Under 
these stabihty conditions, we establish convergence of the hybrid scheme. While 
we only consider the flat interface here, the analysis can be extended to smooth 
interface between the atomistic and continuum regions. In that case, we need to 
check the interface stability condition Assumption [C] for interface with different 
angles. 

For the example of atomistic-to-continuum coupling for triangular lattice consid- 
ered here, the stability can be checked by hand. For many other more complicated 
interactions, this might not be easy or even possible to do by hand. One possible 
direction is to use numerical or symbolic computations to check stability conditions, 
in analogy to checking GKS conditions for example as in [32] . This is an interesting 
future research direction. 

The result in this paper does not apply to transition interface between atomistic 
and continuum regions that involves corners. Extension of the result to hybrid 
schemes with transition interface with corners would be very interesting. 
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